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Discharge current and terminal voltage measurements were performed on a gallium 
electromagnetic thruster at discharge currents in the range of 20-54 kA. It was found that 
the arc impedance has a value of 6-7 mQ at peak current. The absence of high-frequency 
oscillations in the terminal voltage trace indicates lack of the " onset” condition often seen in 
MPD arcs, suggesting that a sufficient number of charge carriers are present for current 
conduction. The mass ablated per pulse was not measured experimentally; however the mass 
flow rate was calculated using an ion current assumption and an anode power balance. 
Measurement of arc impedance predicts a temperature of 3.5 eV which from Saha 
equilibrium corresponds to Z = 2.0 - 3.5, and assuming Z = 2 yields an Isp of 3000 s and 
thrust efficiency of 50%. 
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Nomenclature 


geometric constant, (H/m) V Q 

specific heat (J/kg-K) V a 

elementary charge (C) V s 

enthalpy (J/kg) Z 

discharge current (A) Z arc 

Boltzmann constant (J/K) Z L 

inductance (H) a 

ion mass (kg) 77 

Power (W) r| Ga 

inner electrode radius (m) A 

outer electrode radius (m) X 

time (s) v e i 

temperature (K) a 

exhaust velocity (m/s) p 0 


= capacitor charging voltage (V) 

= anode sheath voltage (V) 

= work function (eV) 

= charge number 
= arc impedance (£ 2 ) 

= dynamic impedance (£ 2 ) 

= current carrying capability of plasma 
= thrust efficiency 
= electrical resistivity (£ 2 -m) 

= impact parameter 
= thermal conductivity (W/m-K) 

= ion-electron collision frequency (Hz) 
= conductivity (S/m) 

= permeability of free space (H/m) 


I. Introduction 

T HE gallium electromagnetic (GEM) thruster 1,2 uses a high-current arc discharge to enable operation at high 
exhaust velocities (15-30 km/s). Gallium appears to offer several advantages over propellants currently in use. 
The non-toxicity and liquid temperature range (30 °C - 2204 °C) provide numerous systems advantages, while the 
high atomic mass and low ionization potential (5.99 eV) predict low frozen flow losses. A previous study 3 has 
presented performance estimates and emission spectroscopy measurements for a 20 J/pulse gallium plasma source. 
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In this paper current and voltage data are presented on a GEM thruster powered by a 1.4 kJ pulse forming network 
(PFN). In the present design, a low energy spark plug pre-ionizes propellant from a center gallium anode, and the 
ablation products expand radially towards the inner walls of the outer cathode. When the gallium bridges the 
electrodes, it effectively acts as a switch causing the main discharge to commence, with the Lorentz force 
accelerating the gallium plasma axially to produce thrust. 

A key issue in the operation of the GEM thruster is that of mass injection. Previous electromagnetic accelerators 
have operated in various modes, which can be classified by the injection scheme of the propellant. The detonation 
(snowplow) acceleration mode is established when the thruster barrel is prefilled with gas, and a moving current 
sheet accelerates the propellant. Gas-fed deflagration MPDTs employ a stationary j x B force distribution, with the 
propellant fed from the back, ionized, and then pumped through the current sheet in a relatively efficient 
acceleration process. In most pulsed plasma thrusters (PPTs), the insulator material is ablated and entrained in an 
accelerating current sheet to produce thrust. The GEM thruster differs from prior EMTs in that the central electrode 
is vaporized by the discharge to provide the ionized material, and the arc thus controls the mass flow rate. A 
summary of various thrusters along with their corresponding injection schemes is given in Table 1. Also tabulated is 
the arc impedance, which is calculated by dividing the arc voltage V arc by the discharge current I. 


Table 1. Comparison of Pulsed Coaxial EM Accelerator Injection Modes 


Accelerator 

Type 

Propellant 

Injection 

Mode 

Fast 

Valve 

Injector/ 

Switch 

Arc 

Impedance, 

mQ 

Marshall Gun 4 

Deuterium 

Gas-puff 

Yes 

No 

6 

Deflagration Gun 5 

Hydrogen 

Gas-fed 

Yes 

No 

50 

Solid-Fed MPD 6 

Teflon 

Ablation 

No 

Yes 

9 

Self-Field MPD 7 

Argon 

Gas-fed 

Yes 

Yes 

7 

Gas-Fed PPT 8 

Argon 

Gas-fed 

Yes 

Yes 

10 

Solid-Fed PPT 9 

Teflon 

Ablation 

No 

Yes 

14 

Lorentz Force Accelerator 10 

Lithium 

Vapor-fed 

No 

No 

59 

Self-Injecting MPD 

Gallium 

Evaporation 

No 

Yes 

6 


Performance estimates for the GEM thruster have been presented in Ref. [3]. Two parameters of particular interest 
are the ion Larmor radius R L and the charge carrier parameter a, which represents the fraction of the discharge 
current that can be carried by the number of ionized particles present 


Ze dm 
ml d t 


( 1 ) 


By examining prior MPD performance data, one can conclude that efficient operation of EMTs occurs only when 
there are enough ion charge carriers present to conduct the current and when the Larmor radius is smaller than the 
thruster channel. For heavier propellants such as gallium this implies: 

• The discharge current must scale with the ion mass and exhaust velocity to keep the Larmor radius 
small. 

• The propellant mass flow rate must scale with ion mass and discharge current to maintain a high 
number of charge carriers, i.e. a > 1. That is, thrusters using heavier propellants such as gallium need 
to be operated at much higher mass flow rates than argon or hydrogen to provide enough current- 
carrying particles. 

Consequently, the present experiments were conducted at relatively high discharge currents (20-54 kA). It was 
assumed that the arc would vaporize the proper amount of gallium to sustain the discharge; that is, the mass flow 
rate was not actively controlled during the experiments. 
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II. Experimental Apparatus 



A. Testing Facilities 

Tests were conducted at the University of Illinois Electric 
Propulsion Laboratory. The thruster was mounted into a 5- 
way vacuum cross, which is attached to the main vacuum 
chamber. A baffle was inserted into the 5 -way cross to 
condense ejected gallium particles before they reach the main 
chamber. Vacuum was maintained by a TPH1500 
turbomolecular pump, which is backed by two Roots blowers 
and two Kinney mechanical displacement pumps. The base 
pressure of the facility during testing was 5x1 0" 5 torr. 

A ten- section, 5 -line pulse forming network (PFN) was 
used to power the thruster. Each section consists of a 88 pf 
capacitors charged to a maximum of 800 V, in series with a 
85 nH inductor, for a total of 50 capacitors yielding a bank 
energy of 1.4 kJ. A copper sheet with punched holes was 
laid over the capacitors to provide a common ground (ground 
plane). Brass standoffs with a height of 3.8 cm were 
machined to provide spacing between the stage inductors and 
the ground plane. The section inductors were composed of 9 
cm long, 0.32 cm (1/8”) diameter copper rod, yielding the 
section inductance of ~ 85 nH. The five PFN lines were 
connected in parallel using a bus made of 2.5 cm wide 
copper strip. To minimize stray inductances in the charging 
circuit, the PFN was mounted directly near the vacuum 

chamber, as shown in Fig. 1. A 7.0 cm wide strip of copper sheet formed the transmission line. The capacitors were 
charged using a Hipotronics DC power supply. 


Figure 1. 3.2 kJ, 31 mQ pulse forming network. 


B. Discharge Initiation 

Arc initiation in a vacuum environment is difficult due to the extremely low vapor pressure of the gallium 
propellant. Prior experiments 3 used a thin layer of carbon to bridge the cathode and anode, creating a conductive 
path for the discharge current. While successful in initiating breakdown, a fraction of the discharge energy was lost 
in excited and ionized carbon atomic states. In the present experiment, a spark plug was used to initiate the 
discharge. The igniter plug used is a semiconductor-type aircraft turbine combustor spark plug mounted 5 cm from 
the face of the thruster. It is powered by a 20 J non-reversing current pulse. It was found that after ~50 shots the 
spark plug became clogged with gallium and would no longer fire. This was remedied by firing the igniter plug 
several hundred times in between tests to keep it clean. 


C. Operating Polarity 

Unlike prior MPDTs, the GEM thruster operates with a central anode. This is in part due to the discharge being 
initiated with a spark plug. Prior PPT studies 9 experienced difficulties igniting the discharge with a central cathode. 
This was attributed to electrons being unable to reach the center electrode, as they were lost to the positive outer 
electrode and/or repelled by the negative center electrode (illustrated in Fig. 2). 
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Figure 2. Left: An insufficient number of electrons are able to reach the center electrode 
when it is negatively biased. Right: Electrons can reach a positively biased center electrode. 

Additionally, the influence of operating polarity on ion trajectory motion in a coaxial, deflagration mode MPDT has 
been discussed in a previous paper. 2 In that work, two cases were presented: a central cathode (the traditional 
configuration for MPDTs) and a central anode (such as is found in PPTs). In the case of the central cathode, ions 
created near the outer electrode enter a region of strong 5-field as they are accelerated radially inward by the electric 
field, and they are deflected axially by the quxB force. For a central anode, ions are created along the axis in the 
region of strong B- field. These ions are more easily deflected axially because they possess a smaller radial velocity, 
and consequently a smaller Larmor radius, than in the case of a central cathode. 

D. Electrode Assembly 

The coaxial gallium pulsed plasma source is shown in Fig. 3. The stainless steel outer electrode has outer and 
inner diameters of 5.72 cm and 2.22 cm, respectively. A threaded Teflon tube fitted inside a boron nitride (BN) 
insulator electrically separates the inner and outer electrode. A 5.0 mm diameter hole drilled through the center of 
the Teflon allows the gallium to make electrical contact with a brass standoff, which in turn is connected to the PFN 
outside of the chamber. Teflon was used to contain the gallium as the boron nitride was too brittle to withstand the 
high pressures at the face of the thruster during operation. 



SS316L Ga 


BN Brass Teflon 



(a) (b) 

Figure 3. a) Photograph and b) cross-section of the coaxial electrode assembly. 
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The inner diameter of the inner electrode was chosen to satisfy two conditions. It has experimentally been found 
that when the outer to inner electrode ratio exceeds 5, MPDTs become prone to the spoke instability. 11 Additionally, 
for proper operation, it is necessary to heat the gallium surface from room temperature to boiling (2477 K) in a time 
comparable to the current pulse rise time. For a wall exposed to a power flux P , for a time r, the temperature 
increase is given by 12 


AT = 


1 


2 {n p c 


1/2 


el 


n 


1/2 

T 


where p , c, A are the density, specific heat, and thermal conductivity of gallium. The power flux P ei is the power 
received by the anode surface due to the incoming flux of electrons, which is composed of 1) the power associated 
with the anode voltage drop V a , 2) the heating due to electron absorption on the anode surface, and 3) the power that 
accounts for the transport of electron thermal energy from the plasma 13 

^ kT 

P el =I{ ys + ~) + \ IV a 

The surface area over which the discharge current acts is currently an unknown. The possible presence of anode 
spots, the movement of such spots, and the convection of heat between the spots complicate this problem. For the 
purpose of this calculation, it is assumed that the current uniformly heats the entire anode surface. This assumption 
yields heating times in the range of 5-10 ps for a 50 kA pulse. 

E. Current and Voltage Measurements 

Current and voltage measurements across the discharge are important as they provide information on the thruster 
impedance, which is related to the efficiency of the thrust generating process. Measurements were obtained for 
capacitor charging voltages in the range of V Q = 300-800 V. A Pearson 4418 current monitor and Tektronix 
P6015A voltage probe were used to measure the discharge current and terminal voltage, respectively. 
Representative traces are shown in Figs. 5 and 6. Note that the PFN impedance is well-matched to the arc 
impedance, minimizing current reversal. 



Figure 5. Discharge current waveforms for a PFN charging voltage of 400, 600, and 800 V. 
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Figure 6. Terminal voltage traces for a PFN charging voltage of 400, 600, and 800 V. 

Several things can be noted from these figures. It can be seen that the “flat top” portion of the current pulse is not 
actually flat; this is attributed to the skin depth effect in the 3.2 mm inductors, and can be resolved in the future by 
using section inductors with a larger cross sectional area. Additionally, there is a lack of high frequency oscillation 
in the voltage trace indicating that the arc is operating below onset and that there are a sufficient number of gallium 
ions present for adequate current conduction at these discharge current levels (a > 1). 


III. Mass Consumption Estimate 

The mass flow rate is calculated using two separate methods. Rearranging Eq. (1), the mass flow rate can be 
calculated directly from the discharge current (assuming a = 1) 


dm ml 
d t Ze 


(4) 


The amount of gallium vaporized and ionized per pulse can also be estimated by considering a power balance at the 
anode surface 14 


P 4- P = P +P + P 

res el vap cond rad 


(5) 


Anode heating takes place by resistive heating P res and electron bombardment P d . The power is dissipated through 
vaporization P vap , conduction P cond and radiation P rad . The problem of a source of uniform current density acting 
over a circular area on the surface of a semi-infinite solid of conduction material has been solved and the resulting 
power input due to Joule heating has the form 15 


0.327' 


VGa 


71 r 


( 6 ) 


where rj Ga is the electrical resistivity of gallium. The power flux P d is given by Eq. (3). The power flux P vap is the 
power dissipated through evaporation of the anode surface 
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( 7 ) 


p =h — 

± vap "tot ^ 

The total enthalpy h tot is comprised of the enthalpy needed for a phase change (enthalpy of fusion + enthalpy of 
vaporization) and then raising the gas temperature to the arc temperature as given by 16 

Kt =hf+h v +c p AT (8) 

To calculate the power dissipated by thermal conduction, it is assumed that the anode material is a semi-infinite 
solid and that the steady state is quickly reached for thermal flow 17 

P cond =n i ' 1 l{T s -T 0 )r i (9) 

where T 0 is the initial temperature (298 K) and X is the thermal conductivity. The power dissipated through 
radiation can be estimated as 


Prad=^ b T>r i 2 ( 10 ) 

where c> sb is the Stefan-Boltzmann constant, which has a value of 5.67 W/m 2 -K 4 . Using equations (5)-(10), one can 
solve for the mass flow rate 


dm 

dt 


0.32/ 2 /7 Ga tnr, +(F a + K + %JTeV_ fc -T 0 )r, - 

hf +h v +c p AT 


rp4 2 
°sb T S * r i 


( 11 ) 


The computed mass flow rates are given as a function of discharge current in Fig. 10. A sheath voltage and electron 
temperature of 6 V and 3.5 eV, respectively, are used. It can be seen that the mass flow rate calculated from an 
anode power balance gives a value slightly less than a doubly-ionized (Z= 2) plasma. An experimentally based 
estimate of Z is given below in Section IV B. Based on these estimates, we adopt Z = 2 pending further 
measurements including spectral data. 



Current, kA 

Figure 10. Comparison of calculated mass flow rates. 
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IV. Discussion 


A. Efficiency Estimate 

The current and voltage data acquired during testing are summarized in Table 2. Also shown in the table are the arc 
impedance and the dynamic impedance Z L , , which is found from 


1 

Z L =-L'u 
4 e 


The relation for the inductance gradient L’ can be found from the EM thrust relation 


18 


T = 






= bl 


2 



where the geometric constant b has a numeric value of 2.02 x 10" 7 H/m for our configuration. 


( 12 ) 


(13) 


Table 2. Plasma Electrical Characteristics at t = 45 jlxs 


PFN 

Charging 
Voltage, V 

Peak 

Current, kA 

Terminal 
Voltage, V 

Arc 

Impedance, 

mQ 

Dynamic 
Impedance, 
mQ (Z= 1) 

Dynamic 
Impedance, 
mQ (Z= 2) 

300 

19.6 

136 

6.9 

0.5 

1.1 

400 

26.0 

180 

6.9 

0.7 

1.5 

500 

34.8 

220 

6.3 

1.0 

2.0 

600 

41.5 

255 

6.1 

1.2 

2.3 

700 

49.5 

303 

6.1 

1.4 

2.8 

800 

54.0 

330 

6.1 

1.5 

3.0 


The exhaust velocity is found by dividing the thrust by the mass flow rate 


T 


u 


e ~ 


m 


(14) 


The thrust efficiency is calculated from the ratio of the dynamic impedance to the arc impedance, assuming that 
thermal heating does not add to the thrust 





L'TIm 

4 Z nvn 


(15) 


The exhaust velocity as a function of discharge current is plotted in Fig. 7, yielding a specific 
discharge current of 54 kA. The efficiency as a function of the specific impulse is shown in 
that the efficiency is a linear function of the specific impulse. 


impulse of 3000 s at a 
Fig. 8. It can be seen 
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Figure 7. Exhaust velocity as a function of the discharge current for an estimated value of Z = 2. 



Isp, s 

Figure 8. Thrust efficiency as a function of the discharge current for an estimated value of Z = 2.. 


In the next subsection, a simple model is used to estimate the electron temperature to support the claim that the 
plasma is doubly ionized. 

B. Electron Temperature Calculation 

Using Gauss’s Law and assuming a spherical geometry, the electric field has a 1/r 2 dependence 

Ci 

E r =^ (16) 

r 


Noting the electric field is the gradient of the electric potential, the integration constant C 7 can be found 


> 

II 

1 

>\£* 

>2 

u 

=>C X = AV 

r ±__L 

r i 


U 


The current density is also assumed to have a 1/r 2 dependence 


( 17 ) 
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J = - 


(18) 


n r 


Ohm’s Law directly relates the current density to the electric field 


j = a 0 (E + uxB) + ^(jxB) 
B 


(19) 


The estimated hall parameter Q is small so the last term is neglected, and assuming u is radial, equations (16) - (19) 
can be combined to yield an expression for the electrical conductivity 


<j n = - 


1 1 


n AVI 


( 20 ) 


The electron temperature is related to the scalar conductivity through the relation 


a 


o 


2 T 312 

— « 1.3xl0 4 — (S/m) 

m e v ei Z In A 


( 21 ) 


where v ei is the electron-ion collision frequency, m e is the electron mass, In A is estimated to be ~ 3, and the electron 
temperature is in units of eV. The conductivity and electron temperature are calculated using the experimental data 
and are shown in Table 3. The potential difference AV and current / used to calculate the conductivity in Eq. (20) 
are taken to be the terminal voltage and peak current values located in Table 2. 


Table 3. Calculated Conductivity and Temperature 


PFN 

Charging 
Voltage, V 

Arc 

Current, 

kA 

Arc 

Impedance, 

mQ 

Electrical 

Conductivity, 

S/m 

Electron 
Temperature, 
eV(Z= 2) 

300 

19.6 

6.9 

14217 

3.5 

400 

26.0 

6.9 

14249 

3.5 

500 

34.8 

6.3 

15582 

3.7 

600 

41.5 

6.1 

16054 

3.8 

700 

49.5 

6.1 

16116 

3.8 

800 

54.0 

6.1 

16142 

3.8 


This data provides a rough estimate of the degree of ionization. Cassibry 19 has calculated equilibrium conditions for 
gallium plasmas and has found Z = 3.5 at 5 Pa and 3 eV. The ionization fraction as a function of pressure is shown 
below in Fig. 9. Calculations at 50 atm yield a Z = 2.0 for T = 3.5 eV. Because the pressure will decrease 
significantly away from the anode, Z is expected to vary from 2.0 to 3.5 in the discharge. We therefore have 
adopted a conservative value of Z = 2 for our preliminary calculations. 
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Figure 9. Gallium ionization fraction as a function of pressure for an electron temperature of 3.5 eV. 


V. Summary 

The results obtained can be summarized as follows: 

• Experiments were conducted at discharge currents in the range of 19-54 kA. At these high currents and 
short pulse lengths, fast valves are not practical, and a method of mass injection such as evaporation 
presents itself as a viable alternative. 

• Arc impedances in the range of 6-7 mQ were measured from the current and terminal voltage. These 
low values compare favorably with prior EMTs and are required for high efficiency. The arc impedance 
slightly decreases with increasing current, suggesting the electron temperature increases with the 
discharge current. Electrical conductivity calculations yield an electron temperature of - 3.5 eV, and the 
degree of ionization is conservatively estimated to be Z = 2. 

• High-frequency oscillations in the terminal voltage trace (“onset”) were absent for discharge currents as 
high as 54 kA, leading to the conclusion that the arc operates below onset for all conditions tested. 

• The mass ablated per pulse was not measured experimentally; however the mass flow rate was 
calculated using two separate methods, both of which roughly agreed. These calculations also suggest 
that ions carry the bulk of the current. 

• The specific impulse and efficiency both linearly scale with the discharge current. Assuming a Z = 2 
plasma results in a calculated Isp of 3000 s and a thrust efficiency of rj = 50%. 
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